clear all; close all; clc;
% Внутренние параметры
R1 = 220; % Ом
C1 = 33e-12; % Ф
L = 100e-6; % Гн
% Ось времени
fmax=1e6;
T = 1 / (1000*fmax);
t = 0:T:(100* 1/fmax);
lt = length(t);
% Выделение памяти
Uc = nan(1, lt);
Ul = nan(1, lt);
dUc = nan(1, lt);
i = nan(1, lt);
i1 = nan(1, lt);
i2 = nan(1, lt);
di1 = nan(1, lt);
% Внешние параметры, внешние фазовые переменные
E = ones(1,lt);
% Начальное состояние
Uc(1) = 0; % Начальное напряжение конденсатора
dUc(1) = 0;
i(1) = 0; % В начальный момент ток в цепи отсутствует
i1(1) = 0;
i2(1) = 0;
di1(1) = 0;
Ul(1) = 0; % Напряжение на катушке
for k = 2:lt
Ul(k) = E(k) - i(k-1)*R1; % Напряжение на катушке
di1(k) = Ul(k) / (L); % вызывает приращение тока,
i1(k) = i1(k-1) + di1(k-1)*T;
i2(k)=( E(k) - i(k-1)*R1-Uc(k-1))/R1;
dUc(k)= i2(k)/C1;
Uc(k) = Uc(k-1) + dUc(k-1)*T;
i(k)=i2(k)+i1(k);
end
figure(1);
plot(t, [E;Uc*1]);
xlabel('t, ns')
ylabel('A')
legend('E(t)', 'Uc(t)')
grid on